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1 Abstract 

This computational study focuses on understanding and quantifying the effects of external 
acoustical perturbations on droplet combustion. A one-dimensional, axisymmetric representation 
of the essential diffusion and reaction processes occuring in the vicinity of the droplet stagnation 
point is used here in order to isolate the effects of the imposed acoustic disturbance. The simulation 
is performed using a third order accurate, essentially-non-oscillatory (ENO) numerical scheme with 
a full methanol-air reaction mechanism. Consistent with recent microgravity and normal gravity 
combustion experiments, focus is placed on conditions where the droplet is situated at a velocity 
antinode in order for the droplet to experience the greatest effects of fluid mechanical straining of 
flame structures. The effects of imposed sound pressure level and frequency are explored here, and 
conditions leading to maximum burning rates are identified. 

2 Introduction 

This study is part of a larger research program designed to understand and quantify the effects 
of pressure oscillations on the combustion characteristics of condensed fuels under microgravity 
conditions [1]. The issue of microgravity acoustic excitation of flames is especially pertinent to 
understanding the behavior of accidental fires which could occur in spacecraft crew quarters. The 
microgravity diffusion flames which could form are known to be very strongly affected by even small 
convective contributions [2]. Such disturbances could also be introduced by pressure perturbations 
within the crew quarters, e.g., as results from ventilation fans or engine vibrations. In normal grav- 
ity, studies have shown that there can be a significant increase in fundamental transport processes 
from reactive interfaces and surfaces with the imposition of an external acoustical field[3-7]. Under 
reduced gravity, natural convection is no longer a dominant process, and the effect of acoustics is 
likely to be far more pronounced. 

2.1 Background 

A comprehensive review of microgravity droplet combustion research through the early 1990s 
may be found in Law and Faeth[8]. A number of studies of droplet combustion in a microgravity 
environment indicate that the classical “ d 2 law” [9] does describe the behavior of the droplet 
diameter, but that flame radius may not follow the law until a significant amount of time and 
droplet size reduction has elapsed[10, 11]. This observation is confirmed by experiments conducted 
in normal gravity but at very low pressures for small droplets[12] and is found to be associated 
with lower concentrations of fuel vapor adjacent to the flame at ignition than during later periods 
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of combustion when greater fuel vapor accumulation has occurred. A variety of fuels have been 
found to exhibit these characteristics [13, 14]. 

In a gravitational field, experiments have shown that there can be a significant increase in fun- 
damental heat and mass transfer rates from reactive surfaces with the imposition of an external 
acoustical field[3-7]. The experiments of Blaszczyk[5] examined the effect of a stationary acous- 
tic wave imposed on a single burning droplet. Blaszczyk measured up to a 14% increase in the 
combustion rate constant with acoustic excitation (compared with no acoustic excitation) in the 
excitation frequency range of 120-300 Hz and sound pressure level (SPL) range 100-115 dB. The 
experiments of Saito, et al.[4] also studied the effects of acoustic waves on single evaporating and 
burning fuel droplets, but with an examination of the effect of the position of the droplet, which 
was situated at either a pressure node or antinode. They found that when the fuel droplet was 
situated at a pressure node (velocity antinode), there was a two to three- fold increase in evaporative 
or combustion rate constants, but when the droplet was located at a pressure antinode (velocity 
node), there was little appreciable change in the evaporation or combustion rates. 

To our knowledge, the only prior study which has involved acoustical excitation for burning 
droplets in microgravity is the very recent work by Tanabe, et al. [15,16]. This experimental study 
involved the examination of burning droplets situated in an open duct, with imposed mean flow, 
and with a loudspeaker situated just outside of the duct to impose acoustic forcing. Fan-generated 
velocities were of the order 20 cm/sec, with acoustic forcing at frequencies ranging from about 50 
hz to above 900 hz and sound pressure levels up to 140 dB. At an SPL of 125 dB, for example, 
the authors record velocities “in the vicinity of the droplet” which oscillate with an amplitude of 
±10 cm/sec, although it is not clear that a standing acoustic wave can be obtained in their duct 
at all flow conditions examined. These researchers observe significant increases in the burning rate 
constant as SPL is increased, especially when SPL increases from 125 dB to 135 dB, where a 50% 
increase in burning rate occurs. They also observe peaks in the droplet burning rate at acoustic 
frequencies near 400 hz and near 900 hz. These observations are consistent with corresponding 
normal gravity experiments by this group, although under microgravity conditions the increases in 
burning rate with increasing SPL are more dramatic than in normal gravity. 

2.2 Present Study 

In the present configuration, the droplet and surrounding diffusion flame are to be situated 
within an essentially one- dimensional (actually cylindrical) acoustic waveguide where standing 
waves axe generated by a loudspeaker (s) placed at the end(s) of the guide, generating acoustic 
perturbations with varying frequency and amplitude. This configuration is shown schematically in 
Figure 1. For a waveguide length of 1 m, standing waves may be generated at frequencies above 200 
hz, for example. Applied sound pressure levels in the experiments may be limited to magnitudes for 
which the droplet shape remains roughly spherical. In the present computations, however, larger 
sound pressure levels will be explored in order to be able to make at least qualitative comparisons 
with the trends seen in the microgravity experiments of Tanabe, et al.[15, 16]. 

The present computational study focuses on the essential changes in droplet burning rate that 
can occur when the droplet is exposed to acoustic excitation. A fully 3D or at least 2D axisymmetric 
representation of the droplet and flame, accounting for droplet and flame deformation, should 
be able to quantify precisely how such flame stretch results in augmented transport processes 
during acoustic excitation. But even a one-dimensional, axisymmetric representation of the essential 
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diffusion and reaction processes occuring in the vicinity of the droplet stagnation point can serve to 
quantify the effects of the acoustics (and associated acoustic time scales) on the reaction processes 
(with associated reaction, diffusion, and strain field time scales). 

Hence as a preliminary study, the present paper describes the effects of perturbations in 
freestream velocity as experienced in the stagnation region of a methanol fuel droplet (methanol 
is one fuel expected to be examined in upcoming microgravity experiments [1]). In order to make 
at least some qualitative comparisons with the microgravity combustion experiments of Tanabe, et 
al.[15, 16], a mean flow U 0 c i is imposed on the droplet so as to always maintain a positive incoming 
velocity toward the drop surface, with or without acoustic excitation. Figure 2 shows a schematic 
diagram of the flow and boundary conditions in the vicinity of the droplet stagnation point. The 
“driving” velocity U Q 0 is considered to be oscillatory in time, with an amplitude representative of 
conditions at a velocity antinode: 

Uoo ~ t/oo 1 UoGfeStTlLO'h (l) 

Uoq2 i n fact can depend on the location at which the “freestream” driver is imposed in the compu- 
tation, if the imposed frequencies are high enough (or wavelengths of excitation are small enough) 
to effectively reduce the perturbation velocity near the droplet. 

Many studies of flames in the stagnation region of a blunt object or flame consider the driving 
flow to correspond to Hiemenz stagnation point flow [17-19], although in reality this is precisely 
accurate for large Reynolds numbers (based on droplet diameter). In the present case, however, the 
velocity perturbations do not create such large Reynolds numbers. For an applied SPL of 140 dB, 
for example, the maximum incoming velocity perturbations produce peak Reynolds numbers of the 
order 50; if the SPL is reduced to 120 dB, the Reynolds number is of the order 5. Thus in order 
to be more consistent with the magnitudes of imposed oscillatory velocities that are experienced 
by the burning droplets, the Oseen solution [17] for low Reynolds number flow about a spherical 
object is used as the farfield “driver” for the flow in the vicinity of the stagnation point, where 
0-*O. 

3 Problem Description and Methodology 

The flow and reaction processes in the vicinity of the droplet are solved according to governing 
equations for species and energy conservation in an axisymmetric stagnation region with a low 
Mach number simplification [20]: 

ft + JS pft(ft)] r = +P pft„(t?)] r + S + C (2) 
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where the heat flux component, neglecting radiative heat transfer, is given by: 


NS 

Qr = AT r + p £ h k D km {Y k )r (5) 

fc=i 

The imposed “strain rate” e is computed from Oseen flow to be, in the case of temporally varying 
inflow Uoo, dependent on both the spatial variable r and time t: 


e (r, t) = 


Uoo 

2r 




(6) 


where r is measured from the center of the droplet and the droplet radius r 0 is assumed fixed at 
present, consistent with a computational reference frame that lies at the droplet surface. 

Mass production rate terms c are evaluated in (4) in standard fashion, assuming an Arrhenius- 
type reaction rate coefficient dependence on temperature. In these computations all physical prop- 
erties of the various species are assumed to be functions of temperature and are obtained using 
a transport package within CHEMKIN II[2 1] . Formulation of the ratio of specific heats for the 
mixture of gases present in a given computational cell is done so as to avoid thermodynamic incon- 
sistencies in the presence of contact discontinuities, as described in Selerland and Karagozian [22]. 
A full kinetic mechanism for the combustion of methanol [23] is used in the present simulations. 
Finally, at the droplet surface in the present simulations, the Antoine relation is used to represent 
equilibrium evaporation. 

The method of operator splitting via fractional steps, after Yamenko (e.g., as described in 
[24]), is employed to solve the system of equations (2). This procedure enables the system to 
be split into two separate equations, one which only includes the advection-diffusion terms and 
one which only includes the reaction rate source terms. The present method uses a third order 
essentially-non-oscillatory (ENO) [25] spatial integration to solve the advection-diffusion equation. 
ENO methods constitute a class of high accuracy, front capturing numerical schemes for hyperbolic 
systems of conservation laws, based on upwind biased differencing in local characteristic fields. The 
present scheme then employs a stiff ODE solver (VODE[26], an improved version of LSODE[27]) 
accompanied by the CHEMKIN-II package[21] to solve the temporally evolving equations with stiff 


source terms. 

The present numerical methodology was validated by computing the steady state combustion 
behavior of a spherical methanol droplet in stagnant air and comparing results with the numerical 
simulations of Marchese, et al [23]. As shown in Figures 3ab, the distributions of methanol fuel, 
oxygen, products, OH radical, and temperature appear to compare reasonably well with these prior 
computations, which used the same chemical mechanism but with a different numerical scheme 
for the axisymmetric droplet. Comparisons are especially accurate in the vicinity of the droplet 
surface and on the “fuel” side of the diffusion flame. Because the present scheme is designed to be 
able to accomodate acoustic time scales, the farfield steady state condition takes much longer (in 
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computational time) to reach in the present computation, and hence there are some discrepancies 
between the present simulations and the steady state simulations of Marchese, et al [23] far from 
the droplet surface. Nevertheless, the essential features of the present model (transport processes 
near the stagnation point) appear to be the same as in these prior studies. 

4 Results 

The present computations were carried out for several different amplitudes and frequencies of 
acoustic excitation of the droplet. In all conditions it is assumed that the droplet is situated at the 
velocity antinode of a standing wave in order to produce the maximum effect on burning, as seen 
by Saito, et al. [4], Hence temporal and spatial pressure variation here is minimal. There is also a 
mean flow Uooi imposed in the flowfield, consistent with the experiments of Tanabe, et al. [15, 16], 
so that there is continuous positive inflow to the droplet stagnation point. 

The computational domain in this problem extends from the droplet surface (with effective 
droplet radius r a = 1 mm) to a distance 2 cm into the flowfield. In most computations shown 
here the mean velocity f/ooi is of the order 22 cm/sec with a sound pressure level of 140 dB, which 
results in a velocity amplitude U<x >2 of 22 cm/sec. In a few cases, both “unforced” flow (in which 
Uooi = Uoo 2 = 0) and flow that is “forced without oscillation” (i.e., = 0) will also be examined. 

As noted previously, the perturbation velocity actually can be reduced at high frequencies for a fixed 
computational domain, since at frequencies above several hundred hertz, the effective freestream 
velocity seen by the droplet will be reduced from its velocity peak within the waveguide. This 
variation in Uoo will be explored in these computations. 

Average fuel mass flux at the droplet surface (temporally averaged over one period) as a function 
of time is shown in Figures 4 and 5, for the cases where I7ooi mid I/qo 2 are both independent and 
dependent, respectively, on the frequency imposed within the waveguide. Quasi-steady burning 
appears to be reached at about 0.12 sec after ignition in both cases. In the case where the freestream 
velocity components are independent of frequency, as in Figure 4, the quasi-steady burning rate 
increases with frequency, but appears to asymptote to the same values for frequencies above 1000 
hz. For the case where the freestream velocity components can be reduced at the edge of the 
computational domain at high frequencies, as in Figure 5, the quasi-steady burning rate increases 
with frequency until about 200-400 hz, then the burning rate decreases with frequency above 400 
hz. Figure 6 displays more clearly the quasi- steady burning rate (at times exceeding 0.15 sec) as 
a function of imposed frequency for this condition. These observations are in fact consistent with 
the observations of Blaszczyk [5] in normal gravity, and with Tanabe, et al. [15, 16] in microgravity 
experiments, who both observed peaks in the droplet burning rate at acoustic frequencies near 
200-400 hz despite the fact that they used different fuels. The fact that the present computations 
produce peaks at nearly the same acoustic conditions suggests that the peaks may be due to the 
reduction in the effective freestream velocity amplitude seen by the droplet at higher frequencies. 
The present simulations also show a monotonic decrease in burning rate as the imposed frequency 
increases above 400 hz, consistent with the observations of Blaszczyk [5], but in contrast to the 
second increase in binning rate seen by Tanabe, et al. [15, 16] at frequencies above 900 hz. This 
second increase may be due to acoustic resonances which are specific to the open duct geometry used 
in these microgravity experiments [15, 16] which is not represented in the present computational 
studies and also not used in the normal gravity experiments of Blaszczyk [5]. 

The instantaneous mass flux at the droplet surface as a function of time (portion of an acoustic 


6 



cycle) for different imposed frequencies is shown in Figure 7a; Figure 7b shows the temporally 
varying freestream velocity Uoo which is imposed for a SPL of 140 dB. Figure 7a represents the 
case where the amplitude of oscillation is considered constant or independent of the computational 
domain. These results show that the excursions in mass blowing/burning rate at the surface are 
large at lower frequencies and diminish at higher frequencies. The mean mass flux for the lower 
frequencies, however, appears to be lower than that for higher frequencies. At higher frequencies 
in fact the mean mass flux at the droplet surface appears to be close to that for the case where the 
incoming velocity is non-oscillatory, i.e., ?/oo2 = 0. This observation is consistent with diffusion as 
well as premixed flame studies that examine the exposure of the flame to an oscillatory strain field 
[22,28]. At very high frequencies, the flame often has difficulty keeping up with the oscillations 
(some of the reactions have characteristic times that become close to the order of the acoustic time) 
and effectively begins to behave like a flame exposed to a constant rate of strain. It is interesting 
to note that there appears to be a phase shift between the temporal mass blowing rate in Figure 
7a and the imposed oscillatory velocity field in Figure 7b; this phase shift increases with higher 
imposed frequencies of oscillation. 

These trends in instantaneous burning rate may also be explained by examining the variation 
in flame (maximum) temperature as a function of the acoustic cycle, as well as the variable position 
of the flame in time, as shown in Figures 8ab. As seen in Figure 8a, the excursions in magnitude 
of the (peak) flame temperature are much greater for lower frequencies, where the oscillations in 
peak temperature at very high frequency are very small and close to the steady state temperature 
computed for a constant inflow condition ( Uoo2 — 0). In many of the lower frequency cases, the 
mean fl am e temperature also appears to be slightly lower than for the steady state case, but this 
is not an entirely consistent trend. The position of the flame shown in Figure 8b explains some of 
these features. At lower frequencies of excitation, the flame periodically moves very close to the 
droplet surface, closer than is achieved at higher frequencies, upon which the mass flow rate peaks 
as seen in Figure 7a. Although the flame temperature itself appears to be at a temporal minimum 
when the flame moves close to the droplet surface, it nevertheless heats the droplet to the extent 
that the evaporative blowing rate increases. The almost indiscernible motion of the flame at very 
high frequencies of excitation (1000 hz and above) confirms that the flame is unable to respond to 
any significant degree to the imposed acoustic field at such conditions. 

Finally, Figure 9 shows the time- averaged fuel mass flux at the droplet surface when the fre- 
quency of the oscillation is fixed at 200 hz but where the amplitude of acoustic oscillation is altered 
such that SPLs vary from 110 dB to 140 dB. These results show that the increases in burning 
rate are not very pronounced as one increases the SPL from 110 to 120 dB, but as one goes from 
120 to 130 to 140 dB (where corresponding velocity perturbations are much higher, according to 
the logarithmic dB scale used for SPL), a significant (25%) increase in burning rate occurs. It 
appears that these large increases in velocity perturbation (at a velocity antinode) can lead to 
especially large increases in the oscillatory rates of strain experienced by diffusion flames bounding 
the droplet, leading to large increases in the droplet burning rate. These observations are consis- 
tent with the experimental observations of Tanabe, et al. [15, 16] in microgravity, although the 
experiments suggested that even higher burning rate augmentation levels are possible. 
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5 Conclusions 


This computational study has focused on understanding and quantifying the effects of external 
acoustical perturbations on droplet combustion. A one-dimensional, axisymmetric representation 
of the essential diffusion and reaction processes occuring in the vicinity of the droplet stagnation 
point is used here in order to isolate the effects of the imposed acoustic disturbance. Thus while the 
present model clearly does not capture all complexities in this acoustically driven, reactive flowfield, 
it does represent the essential phenomena to the extent that results are consistent with experimental 
observations. With increasing imposed frequencies, an increase in the average droplet burning rate is 
observed until a finite range (200-400 hz), after which the burning rate decreases. This phenomenon 
appears to be associated with the effective reduction in the perturbation velocity “seen” by the 
droplet at higher frequencies, and is quantitatively similar to experimental observations in normal 
gravity [5] as well as microgravity [15, 16]. Burning rate is also seen to increase significantly when 
one increases the SPL from 130 to 140 dB, which produces large increases in the oscillatory strain 
rate of the flow, again consistent with experimental observations [15, 16]. 

Clearly, interactions among chemical, diffusion, acoustic, strain field time scales control the 
phenomena observed in this problem. When acoustic time scales are relatively small (at very high 
imposed frequencies), they approach the order of magnitude of some of the chemical times associated 
with combustion reactions, and the reaction itself has difficulty keeping up with the imposed high 
frequency oscillation. At lower imposed frequencies (higher acoustic times), the reaction becomes 
comparatively “fast”, allowing the diffusion-limited flame structure to be perturbed more strongly 
by the imposed acoustic field. 

While future work in this study will include more complex droplet and flame geometries for the 
acoustically driven, burning droplet, it is useful to know that the esstential features of the transport 
phenomena can be extracted from a relatively simple one-dimensional formulation of the problem. 
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8 Nomenclature 

NS number of species 

Q r r-component of heat flux vector 

T temperature 

Y{ mass fraction of species i 

e energy of the mixture per unit volume 

hi enthalpy per unit volume of species i 

p pressure 

r radial coordinate measured from the center of the droplet 
r 0 droplet radius 
t time 

Uqo “farfield” velocity approaching the droplet surface = E/qoI + U^sinut 
v r local gas velocity in the r-direction 

e effective strain rate imposed by farfield Oseen flow 
A thermal conductivity of the mixture 
p mixture density 

to frequency of imposed acoustic oscillation 
toi mass production rate of species i 
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Figure 1: Schematic diagram of a fuel droplet situated at a pressure node (velocity antinode) within 
an acoustic waveguide. 


diffusion 

flame 



Figure 2: Schematic diagram of freestream flow over a fuel droplet and the flow within the stagnation 
region of a burning droplet. 
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Figure 3: Results from the present computational method for steady state species and temperature 
distributions measured from the surface of a spherical methanol droplet burning in quiescent air. 
Comparisons are made with the simulations of Marchese, et al [23]. 
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(a) 



(b) 

Figure 4. Results from the present computational method for the mass flux at the droplet surface 
(burning rate), averaged in time over an acoustic cycle, as a function of time for different imposed 
frequencies. Part (b) displays the results on an expanded scale. Here the freest ream velocity 
perturbation is assumed to be constant, i.e., independent of frequency. The imposed velocity 
perturbation corresponds to a sound pressure level SPL of 140 dB, with a mean incoming velocity 
of 22 cm/sec. 
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Figure 5: Results from the present computational method for the mass flux at the droplet surface 
(burning rate), averaged in time over an acoustic cycle, as a function of time for different imposed 
frequencies. Part (b) displays the results on an expanded scale. Here the freestream velocity 
perturbation is assumed to be reduced at the edge of the computational domain when required by 
a high frequency within the waveguide. The imposed velocity perturbation corresponds to a sound 
pressure level SPL of 140 dB, with a mean incoming velocity of 22 cm/sec. 
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Figure 6: Results from the present computational method for the mass flux at the droplet surface 
(burning rate), averaged in time over an acoustic cycle, as a function of imposed frequency. Here the 
freestream velocity perturbation is assumed to be reduced at the edge of the computational domain 
when required by a high frequency within the waveguide. The imposed velocity perturbation 
corresponds to a sound pressure level SPL of 140 dB, with a mean incoming velocity of 22 cm/sec. 
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Figure 7: (a) The instantaneous mass flux at the droplet surface (burning rate) as a function of 
the fraction of the acoustic period for different imposed frequencies, (b) The imposed velocity 
perturbation at the edge of the computational domain corresponding to a sound pressure level SPL 
of 140 dB, with a mean incoming velocity of 22 cm/sec. 


16 




1755 


b / 


1745 

* 1740 
v 

3 1735 
<0 

i 1730 



forced without oscillation 
/ \ f *25 Hz 

J \ f = 50 Hz 

r \ f = 100 Hz 

*'* f *200 Hz 

1 f =400 Hz 

« £-1000 HZ 

m »■•**, £=2000 Hz 




i J* 1 

v 1 ** 1 1 1 1 * ' 1 1 


.XT’ 


* \ \ 0 PT 

» » a. , ; 

\ ■, \ .# * * 

\ \ Q- ET yr - 


^ ;7 • ■< 


> 0 * ■ < 
* \ V 0" a Pi 

\ * ra . ^ : ■ 

\ ■ \ -s ..*/* 

\ D 0 tf i 


ila kL 
-6'^ 


* • 


V*, ; f 

\ *_ 

\ ■« / 
X / 

V 


X'k, / ; 

\ "■ / / 

X * 


0.25 0.5 0.75 1 1.25 1.5 

fraction of period 


(a) instantaneous flame temperature 


0.16 | h m * * * 


0.15 | 

I 

o 0.14 I 



rr 


ar 


without oscillation — ' — 
f-25 Hz — *4 — 
f = 50 Hz - 
£ = 100 Hz — ©— 
£=200 Hz — 

£ = 400 Hz - -O - 
£ = 1000 Hz 
f -2 000 Hz — 


0.75 1 1.25 

fraction of period 


(b) instantaneous peak temperature location 

Figure 8: (a) The instantaneous peak (flame) temperature in the flowfield as a function of the 
fraction of the acoustic period for different imposed frequencies, (b) The location of the peak 
(flame) temperature in the flowfield as a function of the fraction of the acoustic period for different 
imposed frequencies. The imposed velocity perturbation corresponds to a sound pressure level SPL 
of 140 dB, with a mean incoming velocity of 22 cm/sec. 
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Figure 9: Results from the present computational method for the mass flux at the droplet surface 
(burning rate), averaged in time over an acoustic cycle, as a function of imposed sound pressure 
level SPL at a fixed frequency of 200 hz. Here the freestream velocity perturbation is assumed to 
be reduced at the edge of the computational domain when required by a high frequency within the 
waveguide. 
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